Electronic Structure of Paramagnetic V2O3: 
Strongly Correlated Metallic and Mott Insulating Phase 
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LDA+DMFT, the computation scheme merging the local density approxi- 
mation and the dynamical mean-field theory, is employed to calculate spectra 
both below and above the Fermi energy and spin and orbital occupations in the 
correlated paramagnetic metallic and Mott insulating phase of V2O3. The self- 
consistent DMFT equations are solved by quantum Monte Carlo simulations. 
Room temperature calculations provide direct comparison with experiment. 
They show a significant increase of the quasiparticle height in comparison with 
the results at 1160 K. We also obtain new insights into the nature of the Mott- 
Hubbard transition in V2O3. Namely, it is found to be strikingly different 
from that in the one-band Hubbard model due to the orbital degrees of free- 
dom. Furthermore we resolve the puzzle of the unexpectedly small Mott gap 
in Cr-doped V 2 3 . 



I. INTRODUCTION 

The phase transition between a paramagnetic metal 
and a paramagnetic insulator caused by the Coulomb in- 
teraction between the electrons is referred to as Mott- 
Hubbard metal-insulator transition. 1,2 Reliable micro- 
scopic investigations of this many-body phenomenon are 
known to be exceedingly difficult. Indeed, the question 
concerning the nature of this transition poses one of the 
fundamental theoretical problems in condensed matter 
physics. Correlation-induced metal-insulator transitions 
(MIT) of this type are found, for example, in transition 
metal oxides with partially filled bands near the Fermi 
level. In these systems band theory typically predicts 
metallic behavior. The most famous example is V2O3 
doped with Cr; 3 ~ 5 see Fig. 1. While at low temperatures 
V2O3 is an antiferromagnetic insulator (AFI) with mon- 
oclinic crystal symmetry, the high-temperature param- 
agnetic phase has a corundum structure. All transitions 
shown in the phase diagram are of first order. In the case 
of the transitions from the high-temperature, paramag- 
netic phases into the low-temperature antiferromagnetic 
phase this is naturally explained by the fact that the 
transition is accompanied by a change in crystal symme- 
try. By contrast, the MIT in the paramagnetic phase is 
iso-structural; only the ratio of the c/a axes changes dis- 
continuously. This may be taken as an indication for a 
predominantly electronic origin of this transition. 

To explain an MIT induced by electronic correlations 
one can either investigate a simplified electronic many- 
body model to understand, at least, some of the basic fea- 
tures of the MIT, or employ material-specific approaches 
such as the density functional theory in the local den- 
sity approximation (LDA). Concerning the former ap- 
proach, 1,3 the spin S = 1/2, half-filled, single-band Hub- 
bard model 6-8 is certainly the simplest possible model 
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FIG. 1: Phase diagram of V2O3 showing the MIT as a func- 
tion of pressure and of doping with Cr and Ti; data points 
from McWhan et al. 5 



to be investigated. In particular, the existence of an 
MIT in the paramagnetic phase of the half-filled Hubbard 
model had been investigated already in the early work of 
Hubbard. 6 ' 9 However, while the Hubbard I and III ap- 
proximations 6 ' 9 describe the insulating phase rather well, 
they do not describe a Fermi liquid phase on the metal- 
lic side. On the other hand, the Gutzwiller approxima- 
tion provides a picture of the break-down of the Fermi 
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liquid phase as indicated by the collapse of the quasi- 
particle peak and the simultaneous divergence of the ef- 
fective mass at a critical value of the Coulomb interac- 
tion U c (Brinkman-Ricc transition). 10 However, within 
this framework one cannot describe the Hubbard bands 
which are essential both in the strongly correlated metal- 
lic phase below U c and in the insulating phase above U c . 
With these limitations, the details of the MIT in the 
Hubbard model remained unclear, except for the one- 
dimensional case 11 which is very particular since it al- 
ways describes an insulator, i.e., U c = + . 

During the last few years, our understanding of the 
MIT in the one-band Hubbard model has considerably 
improved due to the development of dynamical mean- 
field theory (DMFT). 12 ' 13 Within DMFT the electronic 
lattice problem is mapped onto a self-consistent single- 
impurity Anderson model. 13 This mapping becomes ex- 
act in the limit of infinite coordination number 12 and al- 
lows one to investigate the dynamics of correlated lattice 
electrons non-perturbatively at all interaction strengths. 
This is of essential importance for a problem like the 
MIT which occurs at a Coulomb interaction comparable 
to the electronic band-width. In particular, DMFT pro- 
vides a framework for deriving a coherent picture of the 
electronic spectrum at all energy scales, i.e., of the inco- 
herent features at high energies (Hubbard bands), 9 and 
the coherent quasiparticle behavior at low energies. 7 ' 14 
At T = 0, the transition from the metallic to the in- 
sulating state is signaled by a divergence of the effec- 
tive mass and the collapse of the Fermi liquid quasipar- 
ticle peak at the Fermi energy for Coulomb interaction 
U -> U~. 13 < 15 - 17 DMFT furthermore revealed the coex- 
istence of the metallic and the insulating phase below 
a critical point at temperature T c , such that there is a 
first order phase transition in agreement with the exper- 
imental observation for V2O3. To investigate the MIT in 
V 2 3 , Rozenberg et al. 18 applied DMFT to the one-band 
Hubbard model. The influence of orbital degeneracy was 
studied by means of the two- 19-21 and three-band 20 Hub- 
bard model for the semicircular density of states (DOS) 
of a Bethe lattice. Most recently, a detailed analysis 22 of 
the conductivity change demonstrated that, except for 
a very narrow region directly at the critical point, the 
critical exponents are of the liquid-gas transition type, in 
accordance with a Landau theory for the Mott transition 
within DMFT. 15 - 23 ' 24 

Although the Hubbard model is able to explain cer- 
tain basic features of the Mott-Hubbard MIT in V2O3 
and its phase diagram, it cannot explain the physics of 
that material in any detail. Clearly, a realistic theory 
of V2O3 must take into account the complicated elec- 
tronic structure of this system. In our previous work, 25 
we therefore applied the LDA+DMFT scheme to study 
the MIT in paramagnetic V 2 3 . 26 ' 27 With LDA spec- 
tra calculated for the crystal structure of metallic V2O3 
and insulating (Vo.962Cr . 038)203 as input for the subse- 
quent three-band DMFT(QMC) calculations, we found 
an MIT, or rather a sharp crossover, at U « 5 eV. 



Due to restrictions in computer resources, the QMC cal- 
culations in Ref. 25 were done at T = 1160 K. Sub- 
sequently, extensive QMC simulations at temperatures 
down to T w 300 K were performed to make possible 
a comparison between experiment and theory at exper- 
imentally relevant temperatures. Those computations 
yielded spectra with a quasiparticle peak at the Fermi 
edge considerably stronger than that at 1160 K and were 
in contrast to the existing photocmission measurements. 
This puzzle was finally resolved by recent improvements 
in photoemission spectroscopy (PES) experiments which 
allowed one to perform high-energy, bulk sensitive PES, 
displaying a prominent peak at Ep in essential agreement 
with the LDA+DMFT results. 28 

In this paper, we provide details of our calclations 
reported in Ref. 28, present LDA+DMFT spectra for 
300 K, 700 K and 1160 K both below and above the 
Fermi edge, and compare them to PES and XAS mea- 
surements. Based on calculations of spin and orbital oc- 
cupations, we then discuss the properties of the ground 
state. In particular, the nature of the MIT turns out to 
be rather distinct from that in a one-band model, i.e., we 
find that the effective mass in the a\ g orbitals does not 
diverge at the MIT transition. 

The paper is organized as follows: In Section II the 
LDA band structure, based on the experimental corun- 
dum crystal structure, and the relevance of V-V pairs 
are discussed. Electronic correlations are taken into ac- 
count by DMFT as described in Section III. The result- 
ing LDA+DMFT spectra are presented in Section IV, 
including a discussion of the dependence on temperature 
and the Hund's rule exchange coupling. The pecularitics 
of the MIT in V 2 3 and the differences to the MIT of 
the one-band Hubbard model are worked out in Section 
V. A detailed comparison with the experimental spectra 
follows in Section VI. A summary and outlook is finally 
provided in Section VII. 



II. CRYSTAL AND ELECTRONIC STRUCTURE 

In the paramagnetic metallic phase stoichiometric 
V2O3 crystallizes in the corundum structure, which has 
a trigonal lattice and space group i?3c (D® d ) with lat- 
tice constants a = 4.9515A and c = 14.003A. 29 Vana- 
dium and oxygen atoms occupy the Wyckoff positions 
(12c) and (18e) with internal parameters zy = 0.34630 
and xo — 0.31164, respectively, which deviate markedly 
from the value 1/3 assumed in an ideal hexagonal ar- 
rangement. 29 Within the corundum structure the vana- 
dium atoms are arranged in pairs along the hexagonal 
c-axis, which can be derived from an ideal chain struc- 
ture by introducing vacancies at every third site. 30 The 
oxygen atoms form distorted octahedra around the vana- 
dium sites. While the V-V pair along the hexagonal c- 
axis shares octahedral faces, the octahedra are interlinked 
via edges and corners perpendicular to this axis, 29,30 see 
Fig. 2. 
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V-V pair 
FIG. 2: Crystal structure of V2O3. 



In the Cr-doped paramagnetic insulating phase the 
lattice symmetry is preserved, but the crystal struc- 
ture parameters change slightly. In particular, for 
(Vo.962Cro. 038)203 the lattice constants amount to a — 
4.9985A and c = 13.912A and the positional parame- 
ters are zy — 0.34870 and xo = 0.30745, respectively. 29 
All these changes combine into a distinct displacement 
pattern: As compared to the metallic phase the shared 
octahedral faces between the V-V pair shrink, while those 
octahedral faces pointing to the opposite side along the 
c-axis, i.e., towards the aforementioned vacancies, are 
enlarged. At the same time, the vanadium atoms shift 
away from the shared faces of the V-V pair such that the 
distances within the V-V pair increase upon Cr-doping, 
even though the c-axis lattice constant decreases. The in- 
creased a-axis lattice constant directly leads to enhanced 
vanadium distances within the a6-plane. As a net result, 
all nearest-neighbor vanadium distances are enlarged for 
insulating (Vo.962Cro. 038)203. Hence, we expect a re- 
duction of the bandwidth especially of the ^g-derived 
bands. From a comparison of pure and doped V2O3 
as well as Cr203, Dernier 29 concluded that the metallic 
properties are intimately connected with the vanadium 
hopping within the afc-plane rather than with hopping 
processes between the V-V pairs along the z-axis. 

In a first step LDA band structure calculations 31 ' 32 
were performed, which used the augmented spherical 
wave (ASW) method. 33 - 34 Figs. 3,5,6, and 7 show these 
bandstructures along selected high symmetry lines (Fig. 
4) within the first Brillouin zone of the hexagonal lat- 
tice and the densities of state (DOS) for V2O3 and 
(V .962Cr . 038)203, respectively. 

In total our results are in good agreement with those 
published by Mattheiss. 30 In particular, while the O 2p 
derived bands show up in the range between -9 and -4 
eV, the V 3d dominated states fall, due to the octahedral 
surrounding with oxygen, into two groups of bands: ti g 
and e g . With this separation, the e g bands will be empty 
and the ti g bands partially filled with two electrons per 
V ion. 

Due to the lower trigonal lattice symmetry the ti g 
states are further split into doubly and singly degener- 
ate e£ and a\ g states, see Figs. 5, 7 and the scheme 
Fig. 8. The value of this splitting (« 0.3 eV for the 
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FIG. 3: Electronic bands of V2O3 along selected symmetry 
lines within the first Brillouin zone of the hexagonal lattice, 
Fig. 4. The width of the bars given for each band indicates 
the contribution from the a\ g orbitals. 




FIG. 4: First Brillouin zone of the hexagonal lattice. 



centers of gravity) is much smaller than the t2 g band- 
width (« 2eV). However, as the value of the Coulomb 
interaction parameter U (U > 4eV) is larger than the 
bandwidth, this small trigonal splitting strongly deter- 
mines the orbital ground state of the V-ion obtained from 
LDA+DMFT calculations, as will be shown below. To 
highlight the difference between e g and a\ g states, we 
append to each band at each k point a bar in Fig. 3, 
whose length is a measure for the contribution from the 
aig orbitals to the respective wave function. 

The changes on going to (Vo.962Cro. 038)203 are stated 
easily: In Figs. 6 and 7, we observe a narrowing of the 
t 2g and e g bands of « 0.2 and 0.1 eV, respectively, as 
well as a slight downshift of the centers of gravity of the 
e g bands. However, the insulating band gap expected 
for a calculation with the insulating crystal structure is 
missing. 

As already mentioned a peculiarity of the corundum 
crystal structure are the c-axis V-V pairs. Long ago 
Allen 35 emphasized the importance of the intra-pair in- 
teractions for interpreting spectroscopic properties of 
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FIG. 5: Total and partial densities of states (DOS) of V2O3 
per unit cell. 
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FIG. 6: Electronic bands of (Vo.962Cro. 038)203. 



V2O3 and its solid solution with Cr 2 03. Since the a ig or- 
bitals are directed along the c-axis, these orbitals are the 
ones which mediate a strong hybridization between V-V 
pairs. This hybridization for the V-V pair led Castellani 
et al. 3e to a model (see Fig. 8) where two of the four 
electrons per V-V pair occupy a bonding molecular or- 
bital formed by a\ g orbitals, leaving two electrons (one 
per site) in a partially filled twofold-degenerate e g band. 
That results in a spin-i orbitally degenerate state per 
V ion with complicated orbital and spin ordering pattern 
explaining the unusual properties of the low-temperature 
antiferromagnetic phase. 

Indeed, Figs. 3 and 5 show some splitting of the a\ g 
bands, in particular, between the M and the K point of 
the bandstructure, as one would expect from the forma- 
tion of a chemical bond. But, the situation is far more 
complicated than a simple chemical bonding of the a\ g 
band into a bonding and an anti-bonding band: There 
is some additional spectral weight near Ep (e.g., in the 
vicinity of the T point). Moreover, there is not even a 
low-lying "bonding" a\ g band in parts of the Brillouin 
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FIG. 7: Total and partial densities of states (DOS) of 
(Vo.962Cro. 038)203 per unit cell. 
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FIG. 8: Left and right: splitting of the t\ g orbitals in the 
corundum crystal structure. Middle: Formation of a chemical 
bond for a single V-V pair along the c-axis. 



zone (e.g. between the H and A point). 

The Castellani et al. model 36 was challenged by Park 
et al. 37 Based on the polarization dependence of x-ray 
absorption experiments they came to the conclusion that 
the V 3+ ion is in a spin-1 state. They also demon- 
strated that the orbital ground state of the ion is pre- 
dominantly e g e g with a small admixture of e g ai g config- 
urations. This was later supported by LDA+U calcula- 
tions of Ezhov et al. 3S where a spin-1 ground state with 
a e g e g orbital configuration was obtained. 

Nevertheless, the picture where the strongest hy- 
bridization parameter in V2O3 is the a\ g — a\ g hopping 
within the V-V pair, with all other hybridizations being 
much smaller, is still popular. Many theoretical studies 
of this material start with a (as good as possible) solution 
for the V-V pair and consider the inter-pair interactions 
as a perturbation. 39-41 All these model calculations were 
based on the values of hopping parameters obtained by 
least-square fit of LDA bands to a model Hamiltonian 
with nearest-neighbor hopping. Recently this problem 
was reexamined by Elfimov et al. 42 who found that the 
value of the a\ g — a\ g hopping in the V-V pair is signifi- 
cantly reduced if next-nearest neighbor hoppings (which 
were found to be significant) are taken into account in 
the fit to a model Hamiltonian. Hence, one cannot con- 
sider inter-pair hoppings as a mere perturbation as it was 
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taken for granted for a long time. 



III. INCLUDING ELECTRONIC 
CORRELATIONS VIA DMFT 

The LDA band structure of the previous section 
clearly fails to describe (Vi-xCr^C^. In particular, 
the chromium doped compound ((Vo.962Cr . 038)203 is 
an insulator whereas LDA predicts metallic behavior. 
The reason for this failure is that LDA deals with elec- 
tronic correlations only very rudimcntarily, namely, the 
dependence of the LDA exchange-correlation energy on 
the electron density is given by perturbative or quan- 
tum Monte-Carlo calculations for jellium, 43,44 which is a 
weakly correlated system. To overcome this shortcom- 
ing, we supplement the LDA band structure by the the 
most important Coulomb interaction terms, i.e., the local 
Coulomb repulsion U and the local Hund's rule exchange 
J. The local Coulomb repulsion U gives rise to a genuine 
effect of electronic correlations, the Mott-Hubbard metal 
insulator transition. 13 ' 15 " 17 ' 23 ' 24 ' 45 If the LDA bandwidth 
is considerably larger than the local Coulomb interaction, 
the LDA results are slightly modified but the system re- 
mains a metal. If the LDA bandwidth is much smaller 
than the local Coulomb interaction one has essentially 
the atomic problem where it costs an energy of about 
U to add an electron and the system is an insulator. 
In between, the Mott-Hubbard metal insulator transi- 
tion occurs with V2O3 being on the metallic side whereas 
(Vo.962Cr . 038)203, which has a 0.1-0.2 eV smaller band- 
width, is on the insulating side. 

Interpreting the LDA band structure as a one-particle 
Hamiltonian H^da an d supplementing it with the local 
Coulomb interactions gives rise to the multi-band many- 
body Hamiltonian 46 



LDA value J = 0.93 eV. 48 Unfortunately, such an ab ini- 
tio calculation is not feasible for the Coulomb repulsion 
U since U depends sensitively on screening which leads to 
uncertainties of about 0.5 eV. 49 For our present purposes 
this uncertainty is too large since V2O3 is on the verge 
of a Mott-Hubbard metal-insulator transition, and, thus, 
small changes of U have drastic effects. In particular, 
due to the small differences in the LDA band structure it 
is unlikely that for a U value calculated by constrained 
LDA, V2O3 is metallic whereas (Vo.962Cr . 038)203 is in- 
sulating. Therefore, we adjust U in such a way as to 
make sure that these two systems are metallic and in- 
sulating, respectively. A posteriori, we will compare the 
adjusted value with those calculated by constraint LDA 
calculations and those extracted from the experiment. 

So far, we did not specify H^ DA . In principle, it should 
contain the valence orbitals, i.e., the oxygen 2p orbitals 
and the five vanadium 3d orbitals per atom and, maybe, 
some additional s orbitals. However, for V2O3 we are in 
the fortunate situation that the three t 2g bands at the 
Fermi energy are well separated from the other orbitals, 
see Fig. 7. Therefore, it is possible to restrict ourselves 
to the three t 2g bands at the Fermi energy which are 
made up of the corresponding atomic vanadium 3d or- 
bitals with some admixtures of oxygen p orbitals. In the 
case of three degenerate t 2g orbitals, which is close to our 
situation where bandwidths and centers of gravity of the 
ai g and the doubly-degenerate e g band are very similar, 
the k-integrated Dyson equation simplifies to become an 
integral over the DOS 27 



H = Hl DA +U° 



(V- S aa J) 4i) 



i m^rh aa 



Here, i denotes the lattice site and h ima is the oper- 
ator for the occupation of the t 2g orbital to with spin 
a £ {Til}- The interaction parameters are related by 
V = U — 2 J which is a consequence of orbital rotational 
symmetry. This holds exactly for degenerate orbitals and 
is a good approximation in our case where the t 2g bands 
have similar centers of gravity and bandwidths. As in 
the local spin density approximation (LSD A), the spin- 
flip term of the exchange interaction is not taken into 
account in Eq. (1). The consequences of this approxima- 
tion for states in the vicinity of the Fermi energy do not 
seem to be large as comparative calculations using the 
non-crossing approximation within DMFT show. 47 Fur- 
thermore, a pair hopping term proportional to J is ne- 
glected since it requires that one orbital is entirely empty 
while another is entirely full which is a rare situation in 
the solid state and corresponds to highly excited states. 
For the Hund's rule coupling J we take the constrained 
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Here G m (u), E m (w), and N^(e) 50 are the Green func- 
tion, self energy, and LDA density of states, respectively, 
for the t 2g orbital to. In principle, A^(e) should contain a 
double counting correction, which takes into account that 
part of the local Coulomb interaction already included in 
the LDA. However, this correction results in the same ef- 
fect for all three orbitals and, hence, only translates into 
a simple shift of the chemical potential fi. This makes the 
issue of how to calculate the double counting correction 
irrelevant for the present purposes. The (shifted) /1 has 
to be controlled according to the vanadium valency, i.e., 
in such a way that there are two electrons in the three 
bands at the Fermi energy. 

Within DMFT the k-intcgratcd Dyson equation (2) 
has to be solved self-consistently together with a one-site 
(mean field) problem which is equivalent to an Anderson 
impurity model with hybridization A m (o/) fulfilling 13 

/CO A / l\ 

do/-^. (3) 
-00 u) — ui 

The self-consistent solution of the Anderson impurity 
model given by (3) together with the Dyson equa- 
tion (2) allows for a realistic investigation of materi- 
als with strongly correlated electrons. At small val- 
ues of U this procedure typically yields a spectrum 



6 



with a central quasiparticle resonance at the Fermi en- 
ergy and two incoherent Hubbard side bands, while at 
larger values of U the quasiparticle resonance disappears 
and a metal- insulator transition occurs. 12 This approach 
has been successfully applied to a number of transition 
metal-oxides, 49,51 transition metals, 52 and elemental Pu 
and Ce. 53 For more details and an introduction to the 
LDA+DMFT approach we refer the reader to Ref. 27. 

In the present paper, we solve the multi-band Ander- 
son impurity model by QMC, 54 where by means of the 
Trotter discretization and Hubbard-Stratonovich trans- 
formations the interacting Anderson impurity model is 
mapped to a sum of non-interacting problems, the sum 
being performed by the Monte-Carlo technique. We em- 
ploy a Trotter discretization of At = 0.25 eV -1 unless 
noted otherwise and follow Ref. 55 for the Fourier trans- 
formation between Matsubara frequencies and imaginary 
time r. To obtain the physically relevant spectral func- 
tion we analytically continue the Green function from 
Matsubara frequencies to real frequencies by means of 
the maximum entropy method. 56 The QMC has the ad- 
vantage of being numerically exact while the main dis- 
advantage is that it is restricted to higher temperatures. 
The room temperature calculations of this paper were 
computationally very expensive, requiring about 200000 
hours CPU time on the Hitachi SR8000-F1 at the Leib- 
nitz Rechenzentrum Munich. For the implementation of 
QMC in the context of LDA+DMFT, including flow di- 
agrams, see also Ref. 27. 
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FIG. 9: LDA+DMFT(QMC) spectra for paramagnetic in- 
sulating (PI) (Vo.962Cro. 038)203 and metallic (PM) V2O3 at 
U = 4.5, 5, 5.5 eV, and T = 1160 K. 



IV. LDA+DMFT SPECTRA 

Using the crystal structure of paramagnetic metal- 
lic (PM) V2O3 and paramagnetic insulating (PI) 
(Vo.962Cro. 038)203, respectively, as input, we performed 
LDA+DMFT(QMC) calculations with one a lg and two 
degenerate bands. The results for the spectra of the 
t 2g bands are shown in Fig. 9. At U = 4.5 eV both 
crystal structures lead to spectra showing metallic be- 
havior, with a lower Hubbard band at about — leV, an 
upper Hubbard band at 4 eV and a quasiparticle peak at 
the Fermi edge (OeV). The peak at about leV is split 
from the upper t 2g Hubbard bands due to Hund's rule 
exchange as we will discuss below. 

By contrast, at U = 5.5 eV, both crystal structures 
lead to spectra showing nearly insulating behavior. The 
lower Hubbard band is strongly enhanced whereas at the 
Fermi edge, a pseudo-gap is formed. Above the Fermi 
energy, only small changes of the two-peak structure are 
visible. 

Apparently, qualitatively different spectra for the two 
crystal structures require an intermediate value of U. 
This is indeed observed at U = 5.0 eV: Whereas pure 
V2O3 now shows a small peak at the Fermi edge (a 
residue of the quasiparticle peak obtained at U = 4.5 eV) 
and is therefore metallic, the Cr-doped system exhibits 
a pronounced minimum in the spectrum implying that 



it is nearly insulating. Due to the high temperature of 
T = 0.1 eV w 1160 K of the QMC simulations one only 
observes a smooth crossover between the two phases with 
a metal-like and insulator-like behavior of the respective 
curves instead of a sharp metal-insulator transition as 
would be expected for temperatures below the critical 
point (i.e., for T < 400 K in the experiment). The value 
of the critical interaction of 5.0 eV is in accordance with 
constrained LDA calculations by Solovyev et aL 48 who 
analyzed the charging energy between di- and trivalent 
vanadium ions in an octahedral oxygen environment for 
LaV03, obtaining a U value for the t 2g orbitals which 
is only slightly smaller than 5 eV. Similar U values of 
4 — 5 eV were obtained by fitting spectroscopy data for 
vanadium oxides to model calculations. 57 

Not only the overall t 2g DOS but also the band- 
resolved spectra of the ai g and e g bands provide valuable 
insight. In Fig. 10, the ai g and one of the two degenerate 
e g spectra are shown at U = 5.5 eV. The basic features of 
the spectrum can be understood as follows: We will show 
in the next section that the predominant local configu- 
ration has two spin- aligned electrons in the e g orbitals, 
i.e., a e g e g spin-1 configuration, with some admixture 
of ai g e g spin-1 configurations. Since there are more e g 
than a\ g electrons, let us for a moment disregard the 
CLi g e g configurations. The lower Hubbard band at about 
— 1.5eV indicates the removal of an e g electron from the 
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(E-E p ) (eV) 

FIG. 10: LDA+DMFT(QMC) spectrum for paramagnetic in- 
sulating (V .962Cro.o38)2 3 (PI) for U = 5.5 eV; T = 1160 K. 



predominantly SgSg spin-1 configurations. In an atomic 
picture (which is a reasonable starting point for the insu- 
lating phase) this e^e^ — > transition leads to an energy 
gain of V — J — /j « — 1.5 eV (the approximate position 
of the lower Hubbard band). On the other hand, the up- 
per Hubbard band describes the e^e^ — > e^e^e^ tran- 
sitions. Since the spin-alignment is lost this transition 
costs an energy U + V — fi = U + J+(V — J — fi) « 4.4 
eV, which roughly agrees with the position of the up- 
per Hubbard band. On the other hand, adding an a\ g 
electron costs 2V - ,u w 2.6 eV or 2V - 2 J - n w 0.7 eV, 
depending on whether this electron is spin aligned or not. 
This Hund's rule splitting results in the two-peak struc- 
ture of the upper a\ g Hubbard band in Fig. 10. Since 
there are also a\ g e^ g spin-1 configurations, one has some 
modifications of this picture, in particular, there is also a 
splitting of the upper band, resulting in a small peak 
at U w 0.5 eV due to ax g e^ — > ai^e^e^ transitions. As a 
consequence of the splitting of the upper Hubbard band 
into a two-peak structure, leading to a peak at about 
0.5 eV above the Fermi edge, the gap in the insulating 
phase is very small, much smaller than V « 3 eV which 
would be expected in a one-band Hubbard model. This 
also explains the puzzle in the attempt to model the op- 
tical gap with a one-band Hubbard model 18 : fitting to 
the small experimental gap one is led to an unrealistically 
small Coulomb repulsion of about 1 eV and a bandwidth 
of less than 0.5 eV. 

To study the metal-insulator transition at experimen- 
tally relevant temperatures, we performed calculations 
at T = 700 K and T = 300 K. Since the computational 
effort is proportional to T~ 3 , those low temperature cal- 
culations were computationally very expensive. Fig. 11 
shows the results of our calculations at T — 1160 K, 
T = 700 K, and T = 300 K for metallic V 2 3 and for 
insulating (Vo.962Cro. 038)203. In the metallic phase, the 
incoherent features are hardly affected when the tem- 
perature is changed, whereas the quasiparticle peak be- 
comes sharper and more pronounced at lower tempera- 
tures. This behavior also occurs in the Anderson impu- 



rity model 58 and has its origin in the smoothing of the 
Abrikosov-Suhl quasiparticle resonance at temperatures 
larger than the Kondo temperature. However, in contrast 
to the Anderson impurity model this smoothing occurs 
at considerably lower temperatures which is apparently 
an effect of the DMFT self-consistency cycle. 
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FIG. 11: LDA+DMFT(QMC) spectra for paramagnetic in- 
sulating (Vo.962Cro. 038)203 and metallic V2O3 at U = 5 eV. 

To study the possible effect of a smaller Hund's rule 
coupling J, we performed additional calculations for a 
reduced value of J — 0.7 eV, keeping V nearly con- 
stant. The results in Fig. 12 show that the positions 
of the upper Hubbard bands are significantly shifted to- 
wards lower energies while the spectra below the Fermi 
energy are hardly affected. This suggests that the phys- 
ical properies do not change much. Indeed, we find, for 
example, that the spin-1 state hardly changes when the 
Hund's exchange is reduced. Even at values as low as 
J = 0.5 eV the local moment stays almost maximal, i.e., 
(m 2 z ) = 3.85 at J = 0.5 eV, implying that unrealistically 
small values of J are required for the Castellani et al. m 
picture to hold. 59 



V. CHANGES ACROSS THE MOTT-HUBBARD 
TRANSITION 

A. Local magnetic moment and orbital occupation 

The spin and orbital degrees of freedom play an im- 
portant role in the paramagnetic phase of V2O3 and 
in the changes occurring across the MIT. For exam- 
ple, we find the squared local magnetic moment \rn\j = 
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(E E p ) (eV) FIG. 13: a) Spin, and b) orbital occupation vs. Coulomb 

interaction U for metallic V2O3. 

FIG. 12: Comparison of the LDA+DMFT(QMC) V 2 3 spec- 
tra at two strengths of the exchange interaction: J = 0.93 eV 
(as obtained from constrained LDA), U = 5.0 eV, V — 
3.14eV, and J = 0.7eV, f/ = 4.5eV, V = 3.1eV. 



((Emt^wT — ^ml]) 2 ^) to nave a value of (ml) » 4, unaf- 
fected by the MIT, see Fig. 13. This value corresponds to 
two spin-aligned electrons in the (ai 9 , e gl , e^ 2 ) orbitals 
and therefore to a spin-1 state in the Mott-Hubbard tran- 
sition regime in good agreement with polarization depen- 
dent x-ray absorption measurements of Park et al. It 
also agrees with measurements of the high temperature 
susceptibility which give the value of Heff = 2.66/xs for 
the effective magnetic moment. 60 This is close to the ideal 
5 = 1 value fj, e ff = 2.83/x.b. Note that when U is reduced 
to U < 3, the Hund's rule coupling J needs to be reduced 
as well to avoid an unphysical attractive Coulomb inter- 
action (namely, a Coulomb energy U— 3 J < would oth- 
erwise be gained when a spin-aligned electron is added 
to a singly occupied site) . It is this reduction of J which 
finally leads to a smaller local squared magnetic moment. 
Our results of a spin state which is essentially unaffected 
by the MIT is in stark contrast to results for the one- 
band Hubbard model where m 2 z changes strongly at the 
MIT 13 (in fact, this quantity had even been used as an 
indicator for the MIT). 

The orbital occupation (Fig. 13) obtained by us clearly 
rules out a a\ g singlet since this would correspond to 
n„, =1, n e 7, — 1. Therefore our results contradict 

the model of Castellani et al. 36 who proposed the for- 
mation of an oi g singlet and hence a spin-1/2 state. At 
all [/-values we find predominantly occupied e g orbitals, 
but with a significant admixture of a\ g orbitals (see Fig. 
13). On the basis of an analysis of their linear dichro- 
ism data Park et al. 37 concluded that the ratio of the 
configurations (e^,eg) and (e g ,ai g ) is equal to 1:1 for 
the paramagnetic metallic phase (PM) and 3:2 for the 
paramagnetic insulating phase (PI). This corresponds to 
an electron occupation of the (a\ g , e gl , e^ 2 ) orbitals of 
(0.4,0.8,0.8) for the PI phase and (0.5,0.75,0.75) for the 
PM phase. At T = 1160 K we find for the PI phase 



("insulating" crystal structure and U — 5.0 eV) occu- 
pations of (0.28,0.86,0.86), while for the metallic phase 
("metallic" crystal structure and U = 5.0 eV) we ob- 
tain (0.37,0.815,0.815). While our results give a smaller 
value for the admixture of a lg orbitals (even more so at 
T = 300 K) , the tendency for the decrease of this value at 
the transition to the insulating state is well reproduced. 
Fig. 13 also shows that in the immediate vicinity of the 
Mott transition the orbital occupation has a considerable 
temperature dependence, with even fewer electrons in the 
a,\ g orbitals at lower temperatures. 

Further experimental evidence for a (e g , e g ) configura- 
tion in the ground state of the V +3 ions in V2O3 comes 
from Brown et al. 61 They measured the spatial distri- 
bution of the field-induced magnetization in paramag- 
netic V2O3 by polarized neutron diffraction. Their re- 
sults show that the moment induced on the V atoms is 
almost entirely due to the electrons in the doubly degen- 
erate e g orbitals with only a minor contribution from the 
a,\ g orbital. For the antiferromagnetic insulating phase of 
V2O3, calculations of the electronic structure by LDA+U 
also yielded a spin-1 ground state for the V 3+ ion and a 
(e g ,e g ) orbital configuration. 38 

The origin for the ground state orbital configuration 
discussed above is easily understood from the LDA DOS 
(Fig. 5) where the center of gravity of the a\ g orbital is 
0.3 eV higher in energy than the corresponding value for 
orbitals. This shift together with the asymmetry of 
the DOS leads to an LDA occupation of about 0.55 for 
the dig and 0.72 for each of the e g bands. The occupation 
of e g orbitals is further enhanced in the strongly corre- 
lated metallic regime and, in particular, in the insulating 
phase where the Coulomb interaction value (U > 5eV) 
is significantly larger than the bandwidth (W ~ 2eV). 
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B. Quasiparticle renormalization and spectral 
weight at the Fermi level 

To study the MIT in detail we have calculated the 
quasiparticle weight Z by fitting a third order poly- 
nominal to the imaginary part of the QMC self-energy 
IrnE(iuj n ) at the lowest Matsubara frequencies lo„ which 
gives Z = (1 — dImYi{iuj) / dioj)^ 1 via the slope of the 
polynominal at uj — 0. The resulting quasiparticle weight 
for the aig- and the e^-bands is shown as a function of U 
in Fig. 14. With increasing U, Z first shows a strong de- 
crease for both types of bands. However, in the vicinity of 
the MIT at about U = 5eV the a\ g quasiparticle weight 
remains constant while that of the electrons goes to 
zero. This behavior of the a\ g quasiparticle weight is in 
striking contrast to the behavior at the MIT in the one- 
band Hubbard model where Z — > 0, such that its inverse, 
the effective mass, diverges. Indeed, from Fig. 14 alone 
one might conclude that the MIT occurs only for the e g 
band. On the other hand, the total LDA+DMFT spec- 
trum clearly shows insulating behavior at large Coulomb 
interactions. 
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FIG. 14: Quasiparticle weight Z for the a\ g and the bands 
vs. U, using the crystal structure of metallic V2O3. 

A quantity which measures the spectral weight at the 
Fermi energy and does not dependent on the analytical 
continuation to real frequencies is given by 



(3 



G(t = (3/2) = / A(w) 



(3 exp{-(3/2 lo) 
7r 1 + exp(— (3uo) 

v 

K(u>) 



duj. (4) 



More specifically, —(3/ttG(t = (3/2) measures A(lo) in the 
region given by the kernel K(lu) which is centered around 
the Fermi energy at to = and has a width proportional 



to T = 1/(3. The results in Fig. 15 show that, for 300 
K, the spectral weight at the Fermi energy disappears 
at a critical value of U between 5.1 and 5.2 eV for both 
types of orbitals in the case of metallic V2O3. These val- 
ues of U agree quite well with the position where the e g 
quasiparticle weight is expected to disappear in Fig. 14. 
For the 300 K data for insulating (Vo.962Cr . 038)203, the 
critical U- value is between 4.9 and 5.0 eV. With increas- 
ing temperature, the MIT is a smeared out to become a 
crossover and, at 1160 K, is only signaled by a change of 
curvature slightly above 5 eV. 

This analysis reaffirms the correctness of our choice 
of the value U = 5eV for the simultaneous descrip- 
tion of the metallic and insulating behavior of V2O3 and 
( V . 962 Cr . 038 ) 2 3 , respectively. 




FIG. 15: Spectral weight of the oi g and orbitals at the 
Fermi energy, as estimated by —(3/ttG(t = (3/2), vs. U. 

We still have to address the question why the a\ g quasi- 
particle weight remains constant across the transition, 
i.e., why the effective a\ g mass does not diverge. The 
DMFT Green function is given by Eq. (2) which under 
the assumption of a Fermi liquid-like self-energy 



S m (w) = i?eS m (0) 



dReYi m (u>) 



du 



and Z = (1 — dReY, m {uj) / du)\ u —o) 1 yields 



G m {w) 



de- 



ZN° m {e) 



Z{^-ReL m {0) -e) 



(5) 



(6) 



Hence, an MIT can either occur if the effective mass 
diverges, i.e., Z — > 0, or if the effective chemical po- 
tential /1 — i?eS m (0) moves, due to electronic correla- 
tions, outside the non-interacting LDA DOS such that 
N^(fi - i?e£ m (0)) = 0. In the case of V 2 3 the latter 
happens as is demonstrated by Fig. 16, where ReH m (0) 
has been approximated by its value at the lowest Matsub- 
ara frequency, ReY, m (uJo) . At the MIT, between U = 5.1 
and 5.2 eV, fi — i?eE aig (0) crosses the upper LDA band 
edge while [i — i?eS e ir (0) moves below the lower band 
edge. 
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This explains the pronounced changes of the orbital 
occupation in Fig. 13 and, in particular, why an MIT 
can occur although the a\ g quasiparticle weight does 
not vanish (Fig. 14). This unexpected feature of the 
MIT has important physical consequences: Since at the 
MIT N° (ji - ReE aig (0)) -> 0, the height of the a lg 
quasiparticle peak goes to zero, rather than its width 
which is given by Z. For the band we have both 
N° w (// - i?eE e j (0)) -> and Z -> such that height and 
width simultaneously go to zero. Therefore the quasipar- 
ticle DOS N° lg (ji - ReT, aig (0))/Z does not diverge. Con- 
sequently, physical quantities proportional to this quasi- 
particle DOS like the linear coefficient of the specific heat 
and the local susceptibility do not diverge, at least for the 
dig bands. 
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FIG. 16: Effective chemical potential \i — iZeE(wo) vs. U. The 
upper and lower band edges of the non-interacting LDA DOS 
are shown as solid lines and the entire LDA DOS of V2O3 is 
plotted vertically at the right z-axis. 



VI. COMPARISON WITH EXPERIMENTAL 
SPECTRA 

To be able to compare with experimental photo- 
emission spectra the LDA+DMFT results were multi- 
plied with the Fermi function at the experimental tem- 
perature (T w 180 K) and broadened with a 0.09 eV 
Gaussian to account for the experimental resolution. 28 
The same procedure was used for the comparison with 
x-ray spectroscopy data (with an inverse Fermi function 
at T = 300 K and a broadening of 0.2 eV taken from ex- 
periment). On the experimental side, the PES of Refs. 
28,62 were corrected for the inelastic Shirley-type back- 
ground which also removes the O 2p contribution. All 
experimental and theoretical curves were normalized to 
yield the same area (which is a measure of the occupation 
of the vanadium ti g bands). 

In Fig. 17, LDA+DMFT results for temperatures 
1160K, 700 K, and 300 K are presented. Besides the 
broad, essentially temperature independent peak at 
about — 1.25 eV corresponding to the lower Hubbard 
band, the three curves clearly show the development of a 
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FIG. 17: LDA+DMFT(QMC) results for the metallic phase 
at different temperatures. 

well-defined resonance-like structure just below the Fermi 
energy when the temperature is decreased. The latter 
peak is what remains of the quasiparticle peak after mul- 
tiplication with the Fermi function. At 1160K it is nearly 
equal in height to the lower Hubbard band, and there re- 
mains almost no minimum between these two features. 

In Fig. 18, the LDA+DMFT results at 300 K are com- 
pared with early photoemission spectra by Schramme 62 
and recent high-resolution bulk-sensitive photoemission 
spectra by Mo et al. 2S The strong difference between 
the experimental results is now known to be due to 
the distinct surface sensitivity of the earlier data. In 
fact, the photoemission data by Mo et a/. 28 obtained at 
hv = 700 eV and T = 175 K exhibit, for the first time, 
a pronounced quasiparticle peak. This is in good quali- 
tative agreement with our low temperature calculations. 
However, the experimental quasiparticle peak has more 
spectral weight. The origin for this discrepancy, for a sys- 
tem as close to a Mott transition as V2O3, is presently 
not clear. 




(E-E p ) (eV) 



FIG. 18: Comparison of LDA+DMFT(QMC) results at T = 
300 K with photoemission data by Schramme et al. 62 and Mo 
et al. 28 for metallic V2O3. 

In Fig. 19 we present the corresponding calculations 
for Cr-doped V2O3: There is a lower Hubbard band 
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centered at about -1 eV as in the metallic phase, but 
a quasiparticle peak at the Fermi energy no longer ex- 
ists. It is interesting to note, however, that there remains 
some spectral weight in the vicinity of the Fermi energy. 
Clearly this is not a Fermi liquid effect, but is due to 
highly incoherent states with a large imaginary part of 
the low-frequency self-energy. With decreasing tempera- 

1 1 1 1 1 1 1 1 1 1 — TV — | r 
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FIG. 19: Comparison of LDA and LDA+DMFT(QMC) re- 
sults for insulating (Vo.962Cr . 038)203. 

tures, this incoherent spectral weight is reduced and is ex- 
pected to vanish for T — > 0. Therefore the resistance in- 
creases with decreasing temperature as is to be expected 
for an insulator. For comparison we also show the LDA 
data in Fig. 19. They give a completely different picture: 
Besides a small peak at about — 0.8eV which is roughly 
in the same energy region as the lower Hubbard band 
of the LDA+DMFT calculations, it shows a strong peak 
slightly below the Fermi energy. Clearly, LDA predicts a 
metallic solution, although the input crystal structure is 
that for insulating (Vo.962Cro.03s) 2 3 . 

While the comparison with PES data provides impor- 
tant insight into the physics of V2O3, more than half 
of the theoretical spectrum lies above Ep. For this re- 
gion we compare our results at 1160 K, 700 K, and 300 K 
with O Is X-ray absorption spectra (XAS) for V2O3 at 
300 K by Muller et al. 63 (see Fig. 20). Since in the XAS- 
data the Fermi energy is not precisely determined, the 
data were shifted so that the peaks at 1.1 eV coincide; all 
curves were normalized to the same area. 

The theoretical spectra above Ep are found to be al- 
most independent of temperature. Just above the Fermi 
energy they all show some structure (i.e., a shoulder at 
higher temperatures developing into a small peak at low 
temperatures (300 K)) which is the residue of the quasi- 
particle peak. Furthermore, at 1.1 eV there is a rather 
narrow peak, and at about 4.2 eV a broad peak. The lat- 
ter two structures are parts of the upper Hubbard band 
which is split due to the Hund's rule coupling J. Hence, 
the relative position of those two peaks can be expected 
to depend sensitively on the value of J. A slightly smaller 
value of J will therefore yield an even better agreement 
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FIG. 20: Comparison of LDA+DMFT(QMC) results with X- 
ray absorption data by Muller et al. 63 for metallic V2O3. 

with experiment. 

The absence of any quasiparticle weight near Ep in the 
XAS data is puzzling. This quasiparticle weight is not 
only present in the theoretical spectra above and below 
Ep, but is also seen in the high resolution PES mea- 
surements by Mo et al. 28 below Ep. This calls for addi- 
tional XAS or inverse photoemission spectroscopy exper- 
iments. For comparison with future experiments, we also 
show the theoretical XAS spectra for Cr-doped insulating 
V2O3 in Fig. 21, where our data have been broadened 
with the experimental resolution of Muller et al. 63 
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FIG. 21: LDA-f-DMFT(QMC) spectra for E > E F for insu- 
lating (Vo.962Cr .038)2 3 . 

The comparison between theoretical and experimen- 
tal spectra for metallic V2O3 is summed up in Fig. 22 
where our LDA+DMFT results for 300 K are shown to- 
gether with the experimental PES data by Mo et al. 28 
and XAS data by Muller et al. 63 To document the theo- 
retical improvement achieved by including the electronic 
correlations with the LDA+DMFT technique we also 
show the results of LDA. We note again that, by ad- 
justing the value of U such that the experimentally de- 
termined crystal structures lead to the correct metallic 
and insulating behavior, the spectrum was calculated 
without any further parameter fit. In consideration of 
this fact the agreement of our results with the experi- 
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mental spectra above and below the Fermi energy is re- 
markably good. Although LDA yields the same gross 
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FIG. 22: Comparison of LDA+DMFT(QMC) results with 
PES data by Mo et al. 2S and X-ray absorption data by Miiller 
et al. 63 for the metallic phase above and below Ef 

features, their weight, position and width neither agree 
with LDA+DMFT nor experiment. The interpretation 
of the two large peaks in the upper half of the spectrum 
is also different within LDA and LDA+DMFT. As de- 
noted in the figure, the peaks in LDA are purely from ti g 
(lower peak) and (upper peak) bands whereas they 
are mainly a\ g for the lower and e g for the upper peak 
in LDA+DMFT, with some admixture of the respective 
other band. 

We note that the e a g bands were not taken into ac- 
count in our calculations. Therefore, while the complete 
LDA curve is normalized to an area of 10 (corresponding 
to ten d-electrons) , the experimental and LDA+DMFT 
curves are normalized to an area of 6 (corresponding to 
the six electrons of the ti g bands). We may estimate the 
position of the e g bands in a LDA+DMFT calculation as 
following: Assuming that the intra-<2g Coulomb interac- 
tion V = U — 2 J w 3eV also applies for the interaction 
between spin-aligned ti g - and e^-electrons, and taking the 
difference between the e g - and t2g-band centers of grav- 
ity of roughly 2.5 eV into account, we expect the e^-band 
to be located roughly at 2.5 eV+3 eV= 5.5 eV above the 
lower Hubbard band (-1.5 eV), i.e., at about 4eV. With 
this estimate we expect the (upper) XAS peak at around 
4eV in Fig. 22 to be an admixture of e" g and ti g states. 
More precisely, this upper Hubbard band describes tran- 
sitions from e g e g configurations with two electrons to 
the three-electron configurations e g e g e g , e^e^a lg , and 
e g e g e g (there is a minor admixture of e g a\ g states). 

The properties of paramagnetic V2O3 across the MIT 
obtained with LDA+DMFT for a multi-band model are 
thus found to be remarkably different from those known 
from the one-band Hubbard model. Indeed, the orbital 
degrees of freedom are seen to play an essential role: They 
are not only responsible for the high asymmetry of the 
spectra below and above the Fermi energy, but are also 
required to explain the reduction of the height of the 
quasiparticle peak at the Fermi energy when the MIT is 
approached in the metallic phase, as well as the smallness 



of the insulating gap. 

VII. CONCLUSION 

Using LDA-calculated densities of states for param- 
agnetic metallic V2O3 as well as paramagnetic insulating 
(V .962Cro.o38)2 3 as input, we performed DMFT(QMC) 
calculations at 300 K, 700 K, and 1160K for various U 
values. For U ~ 5eV, the calculated spectra show 
a Mott-Hubbard MIT (or rather a sharp crossover at 
the temperatures accessible by present-day QMC sim- 
ulations). The details of this MIT are quite differ- 
ent from those obtained within the one-band Hubbard 
model. 13,15-18 In the latter model the height of the quasi- 
particle peak at the Fermi energy is fixed and the MIT is 
signaled by a divergence of the effective mass (or the in- 
verse quasiparticle weight 1/Z) such that the width of 
the quasiparticle peak goes to zero. In contrast, our 
LDA+DMFT results show that, for the a lg quasiparticle 
peak, the height goes to zero while the width stays con- 
stant, as indicated by a roughly constant value of 1/Z at 
the MIT. For the e g quasiparticle peak a combination of, 
both, a reduced height and width at the MIT is found. 
This new type of physics, but also the high asymmetry 
of the spectra below and above the Fermi energy as well 
as the smallness of the insulating gap, are all due to the 
orbital degrees of freedom. 

We compared our theoretical data at U = 5eV with 
the results of various experimental measurements and 
found the orbital occupation to be predominantly of e g 
character (with a small admixture of a\ g ) in agreement 
with experiments. The occupation decreases for higher 
L7-values, especially at low temperatures. Furthermore, 
we found a spin-1 state across the MIT in agreement with 
polarization dependent X-ray absorption measurements. 
The 300 K spectrum calculated by us for metallic V2O3 
is in good overall agreement with new bulk-sensitive PES 
measurements. 28 On the other hand, the difference in the 
quasiparticle weight remains to be explained. The com- 
parison with X-ray absorption measurements shows that 
our LDA+DMFT(QMC) calculations also give a good 
description of the spectrum above the Fermi energy. 

All calculations described above were done using the 
integral over the LDA density of states (DOS) (equation 
(2)) to obtain the lattice Green function. For a non-cubic 
system, this procedure is an approximation to the exact 
LDA+DMFT scheme. In the future we plan to make use 
of the full Hamiltonian H° (eq. (1)). In this way it will 
be possible to study the influence of correlation effects 
on all orbitals including the e g orbitals and the oxygen 
states. 
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